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Abstract 

The  traditional  GPS  common  view  technique,  using  C/A  code  receivers,  is  the  main  time 
transfer  method  used  by  various  timing  laboratories  over  the  world.  Moreover,  this  method 
is  used  to  realize  the  TAI  (Temps  Atomique  International)  and  the  TA  (F)  (Temps 
Atomique  Francois).  Clock  offsets  between  laboratory  docks  are  determined  according  to 
a  fixed  procedure  defined  by  the  CCTF  (Comite  Consultatif  du  Temps  et  des  Frequences). 
Using  this  procedure,  one  can  perform  on  average  54  tracks  per  day  (theoretically  90  tracks 
per  day),  providing  dock  offsets.  Each  of  these  dock  offsets  results  from  one  780 ’s  track 
and  is  obtained  as  a  result  of  a  quadratic  regression,  followed  by  various  model-based 
corrections  and  finally  a  linear  regression.  The  dock  offsets  are  then  issued  with  their 
standard  deviations.  However,  this  simplified  estimate  does  not  take  into  account  the 
statistical  properties  of  the  different  types  of  noise  present  in  the  measurement.  We  propose 
here  to  rigorously  estimate  this  time  uncertainty  for  various  types  of  noise  that  characterize 
the  transmitted  GPS  time  offset  data.  This  is  achieved  by  the  calculation  of  the  covariance 
matrix  of  time  samples.  This  method  provides  us  with  the  variances  of  the  drift  coefficients 
and  of  the  residuals,  in  the  case  of  a  linear  drift  model  for  1-day  sample  sets,  taking  into 
account  the  different  types  of  noise.  In  this  paper,  we  compare  the  obtained  results  with 
simulated  and  real  data  over  several  days. 


I.  GPS  TIME  TRANSFER 

The  International  Atomic  Time  (TAI)  scale  is  computed  by  the  Bureau  International  des  Poids  et 
Mesures  (BIPM)  from  a  set  of  atomic  clocks  distributed  in  several  timing  laboratories  around  the 
world.  The  time  transfer  procedure  presently  used  for  the  realization  of  TAI  is  based  on  the  common 
view  approach  [1J.  The  principle  to  compare  remote  clocks  for  the  computation  of  TAI  is  to  connect 
each  timing  laboratory  clock  to  a  GPS  receiver  and  to  simultaneously  observe  the  same  satellite. 
Thanks  to  a  simultaneous  observation  of  the  same  satellite  by  two  remote  time  laboratories,  the  time 
offset  between  each  station  clock  and  the  satellite  clock  is  deduced  by  simple  subtraction  according  to 
the  laboratory  and  satellite  positions. 

About  50  tracks  are  realized  per  day.  To  compute  the  TAI,  the  participating  time  laboratories  have  to 
provide  one  local-to-GPS  time  offset  value  per  day.  This  value  is  deduced  from  a  linear  interpolation 
of  time  offset  samples  accumulated  over  1-day  tracks  (Figure  1).  The  daily  time  offset  if  issued  with 
an  estimation  uncertainty  associated  to  the  interpolation. 
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II.  TIME  DEVIATION  DATA  INTERPOLATION 

The  first  step  when  processing  the  deviation  data  is  to  carry  out  a  linear  interpolation  of  the  N  samples 
{x0,  xx,...,xN_x}  we  have  gotten  every  day  at  the  moments  {t0,  tv...,tN_ J.  This  permits  us  to 
estimate  the  daily  time  offset  and  also  the  uncertainty  of  this  estimate.  The  method  usually  used  to 
obtain  the  interpolated  function  g(t)  =  at  +  b  of  this  samples  is  the  quadratic  least-squares 

interpolation.  The  uncertainty  is  then  obtained  by  estimating  the  daily  standard  deviation  (7 r)  as 

ff2o=  (!) 

W  ^  1=0 

Interpolation  Using  Chebyshev  Polynomials 

Rather  than  this  method,  we  use  the  first  two  Chebyshev  polynomials  [2-4]  as  interpolating  functions. 
The  interpolated  function  x  (t)  of  the  time  deviation  data  is  gotten  as: 

x(t)  =  p{)(p{)  (0  +  A#  (0  +  «(0  (2) 

where  (/)„ (l) ,  <p[  (l)  are  the  first  and  the  second  Chebyshev  polynomial  and  e (l)  is  the  error  function 
of  the  interpolated  function  x (t) .  This  error  function  assumed  random  behavior  of  x ( t ) .  All  the 
parameters  p()  and  px  have  the  same  dimension  as  x(l )  and  are  to  be  estimated  by  p0  and  /), .  As 
shown  below,  this  estimate  is  greatly  simplified  by  orthogonality  and  normality  properties  of  the 
Chebyshev  interpolating  functions. 


Estimation  of  the  Parameters 

Let  us  define  the  vector  <T> ,  associated  with  the  interpolating  functions  (/)  ( I )  as 


O 


f  rt/(^o) A 


yrt/(^Ai- 1)  j 


The  matrix  [O]  is  built 


as 


[o]  = 


rti(^o)  rt(^o) 


vrto(^AT-l)  rt(^AT- l)y 


Let  us  define  now  the  vector  X  as 


X  = 


\XN-\  J 


(3) 


(4) 


(5) 


According  to  (2),  we  may  express  the  vector  X  as 
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X  =  [0\P  +  E 


where  P  is  the  vector  of  the  two  parameters  we  have  to  estimate,  and  E  the  error  vector  standing  for 
the  purely  random  part  of  X: 


Vo" 

and 

E  = 

'  e(t0) 

v  A  / 

l)y 

The  parameter  vector  P  may  be  estimated  using  this  formula 

[<D]TX  =  [0]T[0]P  +  [0]T£  (8 

The  orthonormality  property  of  the  Chebyshev  polynomials  implies: 

[0]T[0]  =  [/3]  (9 

In  addition,  because  of  the  overall  average  (N1  E  —  0^1  the  parameter  vector  may  be  estimated  as: 

P  =  [<D]TX  (10 

Thus,  the  estimate  p .  of  the  parameter  p.  is  obtained  by  calculating 

=  (11 


III.  UNCERTAINTY  ESTIMATION 

The  residuals  may  be  defined  from  ( 1 0)  as  a  vector  8 : 

S  =  W[<S>]P 

The  variance  of  the  residuals  O]  may  be  estimated  by 

(ST-S) 

N'  1 

Knowing  that  the  overall  average  (F[*]p)  =  0  and  using  ( 12),  we  have 

(st-s)=(xt-x)-(pt-p) 


(12) 


(13) 


(14) 


One  of  the  main  properties  of  the  Chebyshev  interpolation  parameters  is  that  cov  ( p0  ,  px )  =  0  for 
any  type  of  noise.  The  variance  of  the  residuals  may  thus  be  estimated  as 
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O, 


_2  .  2 

„+<Ta 


(15) 


where  the  variance  o\  of  the  x{t)  is  equal  to  the  scalar  product  — ^ XT  •  Xj .  As  for  the  scalar 
product  ^/'T  ■  ,  it  is  equal  to  the  sum  of  the  variances  of  each  estimate  p0  and  px . 


Estimation  of  the  Variance  of  the  Residuals  from  the  Noise 
Levels 

Correlation  of  Samples 

We  have  to  estimate  the  uncertainty  of  time  offsets  between  two  time  laboratory  clocks  every  day. 
This  estimate  must  be  evaluated  according  to  the  level  of  the  various  noises  that  are  in  the  time 
samples.  The  autocorrelation  function  R  (l )  of  these  time  samples  x(t )  contains  information  about 

the  type  and  the  noise  levels.  The  power  spectral  density  Sx(f)  is  the  Fourier  transform  of  the 

autocorrelation  function  Rx(t) .  The  two-sided  Sxs(f)  and  the  one-sided  SJf)  power  spectral 
densities  are  thus  defined  as 


JX(/)  =  2-5f(/)  if  f>  0 

k(/)=o  if  /<  o 

so  we  have 

Rx(t]-ti)  =  {x(ti)-x(tj)) 

=  \s?{f)-ei2,cfit‘-,j)-df 

+oo 

=  \Sx{f)-Cos[2xf-(t]-ti)\-df 

0 


(16) 


(17) 


By  taking  into  account  the  low  cut-off  frequency  f)  and  the  high  cut-off  frequency  fh ,  we  may 
rewrite  ( 1 7)  as 


R(!]-ti)=  \ Sx (/) •  cos [2;r/ •  f )] •  df  (18) 

fl 

We  modeled  the  different  types  of  noise  according  to  the  power-law  model  of  Sx(f)  as 

(19) 

a=- 4 


The  different  types  of  noise  are  identified  by  the  values  of  a  as 
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•  white  PM  for  a  =  0, 

•  flicker  PM  for  a  =  - 1 , 

•  white  FM  for  a  =  -2, 

•  flicker  FM  for  a  =  -3, 

•  random-walk  FM  for  a  =  -4. 

Table  1  shows  the  noise  levels-based  formulas  of  the  autocorrelation  function  Rx(t )  and  also  the 
intercorrelation  function  Rx(tj  -  f )  of  time  offset  sequence  [4]. 


Table  1.  Correlations  of  the  time  offset  data  x{t )  versus  the  noise  levels  ka. 


off  frequency. 

Variances  of  the  Parameters  and  the  Residuals 

As  shown  in  (15),  we  have  to  calculate  the  variance  of  the  residuals  from  the  variance  of  the  time 
offset  samples  x  ( t )  and  that  of  the  parameters  p0  and  /), .  Thus,  we  have  to  know  (T\  and  the  <J2p  . 

The  variance  of  x(t)  is  given  as 

={x2(t))  =  Rx(t)  (20) 

Le  us  consider  the  covariance  matrix  of  the  parameters  as 

[Cj  =  (P-PT)  =  [0]T(V-VT)[0]  (21) 


Considering  the  element  whose  position  in  the  matrix  is  (i,  j),  we  may  write 


(x-xT) 

.=Wj-0 

(22) 

thus 

[cX 

'J  k= 0  ‘  1=0 

(23) 
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N-\N-\ 

then  (P.-Pj)  =  [Cv][  =  -  tk)  ■  ■  <p,(tk)  (24) 

k=Q  1=0 

and  so  G2pj  =  (P2 }  =  [Cv ]| .  J=.  (25) 


Finally,  the  variance  of  the  residuals  is  thus  calculated  according  to  (15)  with  the  values  of  the 
variance  of  the  time  offset  samples  given  by  (20)  and  also  with  the  variances  of  the  interpolating 
parameters  (Table  2)  shown  in  (25).  This  uncertainty  estimation  of  the  interpolated  time  offset  value 
is  realized  using  only  the  noise  levels,  since  it  is  based  on  the  autocorrelation  function  of  the  time 
offset  samples. 


Table  2.  Variance  of  the  interpolation  parameters  G 2  and  G 2  versus  the  noise 
levels  ka  [4], 


The  coefficients  {k0,...,  k_4}  represent  the  noise  levels  of  the  power  spectral  density  of  the 
time  offset.  C-0.5772  is  the  Euler  constant.  T0  is  the  time  offset  sampling  period.  The 

low  cut-off  frequency  f)  is  assumed  to  be  much  lower  than  ^  .  Assuming  a  sampling 
satisfying  the  Shannon  rule,  the  high  cut-off  frequency  is  fh  =—  . 

2rn 


Experimental  Validation 


Coherence  Assessment  Using  Simulation 

Aiming  at  realizing  the  assessment  of  the  uncertainty  estimation  given  by  (15),  (20)  and  (25),  we  first 
identified  the  different  noise  types  the  time  offset  samples  actually  contains.  Then  we  generated  a  set 
of  sampled  signals  in  accordance  with  the  levels  of  these  noise  types. 

Figure  2  shows  a  real  data  sequence  taken  as  an  example  to  be  analyzed.  It  represents  time 
connection  data  between  the  Observatoire  de  Besanqon  and  the  Observatoire  de  Paris  between  7 
September  2003  and  2  December  2003.  From  this  signal,  we  generated  1000  time  deviation 
sequences  of  26  days.  We  determined  the  levels  of  the  different  types  of  noise  as  in  [3,6].  According 
to  the  obtained  results,  only  white  FM  and  white  PM  with  the  levels  found  are  included  in  the 
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simulated  sequences. 

To  compare  the  classical  uncertainty  estimation  CD  as  in  (1)  and  the  noise-based  uncertainty  <JE 

given  by  (15),  we  reflected  the  value  of  the  noise-based  standard  deviation  we  calculate  over  26  days 
on  all  of  the  set  of  daily  standard  deviations.  As  shown  in  Figure  3,  the  new  uncertainty  estimation 

<?2e  we  calculated  over  a  fixed  number  of  days  agrees  with  the  daily  uncertainty  values  <JE . 

Figure  4  represents  the  dispersion  of  the  noise-based  standard  deviation  <JE  and  that  of  the  daily 
standard  deviation  <JD.  This  result  is  obtained  from  the  1000  sequences  we  generated.  These 

standard  deviations  show  a  %2  law.  Flowever,  the  dispersion  of  the  noise-based  variance  is  about  5 
times  thinner  than  the  dispersion  of  the  statistic -based  one.  We  may  thus  conclude  that  the 
noise-based  uncertainty  estimation  is  a  better  indicator  than  the  daily  uncertainty. 

Taking  into  consideration  the  % 2  distribution  of  the  uncertainty  dispersion,  we  rate  the  value  interval 
where  the  daily  standard  deviations  should  be  included.  For  example,  in  Figure  4  <JE—  2.68  ns  and 
the  number  of  degrees  of  freedom  1  =  90.  The  confidence  interval  of  90%  is  determined  [5]  to  be 

0.812  •  <(T2d  <1.336  crl  (26) 

According  to  the  value  of  <7E  ,  we  calculate  the  confidence  interval  of  90%  to  be  estimated  as 

2.411  <aD  <3.093  (27) 

The  result  is  represented  in  Figure  5  in  the  case  of  one  example  among  the  1000  sequences  we 
simulated.  It  shows  again  the  good  agreement  between  the  daily  uncertainty  estimation  and  the 
noise-based  one  via  the  deduced  confidence  interval. 

Real  Time  Transfer  Data 

We  operated  the  noise-based  uncertainty  estimation  on  the  real  data  we  used  above  to  determine  the 
typical  noise  levels.  The  results  are  shown  in  Figure  6.  We  notice  that  the  time  transfer  uncertainty 

based  on  the  noise-based  variance  (T2E  is  relatively  underestimated  compared  to  the  daily  uncertainty 
that  is  calculated  as  in  (1).  The  difference  is  evaluated  on  average  to  be  15%. 

Several  differences  between  the  simulated  data  and  the  real  time  offset  samples  may  explain  this 
difference:  On  the  one  hand,  some  samples  are  actually  missing  sporadically,  while  the  simulated 
signals  are  always  complete.  On  the  other  hand,  there  are  some  occasional  high  dispersions  of  time 
offset  evaluation.  These  may  explain  why  the  estimation  of  the  daily  standard  deviation  CD  has  been 
overestimated. 

IV.  CONCLUSIONS 

On  GPS  time  transfer,  the  quadratic  least-squares  interpolation  is  the  usual  method  to  estimate  the 
daily  time  offset  between  two  time  laboratory  clocks.  This  time  offset  information  is  issued  with 
other  significant  information,  namely  the  time  offset  uncertainty.  That  information  is  obtained  by  the 
calculation  of  the  daily  standard  deviation  of  the  residuals  of  the  interpolated  samples. 

In  this  work,  we  chose  to  perfonn  Chebyshev  interpolation  in  order  to  lead  to  a  formulation  of  the 
time  transfer  uncertainty  according  to  noise  levels  contained  in  the  GPS  time  offset  data.  We  saw  that 
the  results  show  a  good  agreement  with  the  simulated  time  offset  signals.  They  also  proved  that  the 
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noise-based  uncertainty  estimation  over  a  few  days  is  a  better  indicator  than  the  daily  uncertainty 
estimation.  However,  we  noticed  that  in  the  case  of  real  time  offset  samples,  the  uncertainty  was 
underestimated  compared  to  the  set  of  daily  standard  deviations.  This  discrepancy  is  probably  due,  on 
the  one  hand,  to  a  few  sporadically  missing  samples  which  are  not  taken  into  account  in  the 
simulations.  On  the  other  hand,  the  connection  discontinuity  and  perhaps  also  the  occasional 
dispersion  imply  an  increase  in  the  daily  standard  deviation  when  the  noise-based  one  is  estimated  in 
accordance  of  the  noise  levels  over  several  days. 

As  the  noise-based  estimation  agreement  with  the  classical  one  is  checked  in  this  work,  that 
discrepancy  does  not  question  the  truthfulness  of  the  noise-based  uncertainty  we  fonnulated.  But  it 
shows  the  necessity  to  evaluate  the  inherent  uncertainty  offset  due  to  the  occasional  dispersion  and 
discontinuity  and  also  the  necessity  to  take  into  account  missing  samples. 
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GLOSSARY  OF  SYMBOLS 


(j  Standard  deviation  of  the  residuals  estimated  statistically  over  1  day 
Number  of  samples  in  a  1  -day  track 
Linear  interpolation  function  of  the  time  offset  samples  x;. 
time  offset  samples  at  the  moments  tt 
Linear  interpolation  parameters  of  g(t) 

Sampling  moments 

Time  offset  sample:  the  time  difference  between  the  remote  and  local  clocks 

with  is  {0,1}  .  The  1st  and  2nd  Chebyshev  polynomials,  i  is  the  degree  of  the  polynomial 
with  i  e  {0,1} .  The  first  and  second  Chebyshev  parameters 
Purely  random  behavior  of  the  time  offset  x(l  ) 

with  i  g  {0,1} .  Estimates  of  the  first  and  second  Chebyshev  parameters 
Vector  whose  N components  are  the  values  of  /(/)  at  each  measurement  time 


N 

g(t) 

x,. 

a,  b 

t, 

x(0 

<PM) 

Pi 

e(t) 

A 

[o] 

p 

E 

X 

W 

P 

S 


Matrix  of  two  column  vectors  O0  and  Oj 

Vector  whose  components  are  the  first  and  the  second  Chebyshev  parameters 

Vector  form  of  e(t) .  The  components  of  E  are  the  values  of  e(t)  for  each  sample  x; 
Vector  form  of  x(t ) .  The  components  of  X  are  the  N  time  offset  samples 
The  unit  matrix  m  x  m 

Vector  whose  components  are  the  estimates  of  the  1st  and  2nd  Chebyshev  parameters 
Vector  whose  components  are  the  Vresiduals  of  the  interpolated  sequence 
Standard  deviation  of  the  residuals  estimates  in  accordance  with  the  noise  levels 


(j  Standard  deviation  of  the  time  offset  sample  x(t) 


cr 


Standard  deviation  of  the  Chebyshev  parameters  /?, 


Sx(f )  One-sided  power  spectral  density  of  the  time  offset  sequence  x(t) 

S?(f)  Two-sided  power  spectral  density  of  the  time  offset  sequence  x(t) 

R  (t )  Autocorrelation  function  of  the  time  offset  sequence  x(/) 
j]  Low  cut-off  frequency  of  the  spectral  density  S  x{f) 
fh  High  cut-off  frequency  of  the  spectral  density  Sx(f ) 

ka  Noise  level  associated  with  the /“noise  of  the  power  spectral  density  of  time  offset  S  x  (f) 
[Cv]  Covariance  matrix  of  the  Chebyshev  parameters  p0  and  px 
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Figure  l.Time  offset  signal  between  two  time  laboratory  clocks.  Dashed  line  represents  the  daily  interpolation  line  usually 
performed  to  estimate  the  daily  time  offset  issued  with  de  the  daily  standard  deviation. 
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Figure  2.  Real  time  offset  samples  between  the  Laboratoire  de  Besangon  and  the  Jaboratoire  de  Paris  (September  seventh  2003  to 
December  second  2003).  This  signal  is  used  as  a  model  to  generate  about  one  thousand  simulated  time  offset  sequences.  They  have 
the  same  noise  levels  as  this  real  signal. 
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Figure  3.  Assessment  of  the  noise-based  estimation  of  time  transfer  uncertainty  for  a  simulated  sequence.  Cross  (+)  represent 
simulated  time  offsets  (eg  over  5  days).  The  boxes  with  its  error  bar  represent  the  daily  interpolated  time  offset  and  its  standard 
deviation.  Dashed  lines  represent  the  noise-based  standard  deviation  around  the  interpolation  solid  line. 
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Figure  4.  Dispersion  for  one  thousand  sequences  of  the  noise-based  standard  deviation  (in  bold  gray  boxes)  and  the  daily  standard 
deviation  (in  hollow  boxes).  Dashed  line  represents  the  theoritical  noise-based  standard  deviation  according  to  the  noise  levels 
entered  for  the  simulation. 


35th  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


05 


(su)  uopeiAap  piEpiiE)s 


149 


-97540 


35th  Annual  Precise  Time  and  Time  Interval  (PTTI)  Meeting 


(SU)  13SJJO  3UI11 


150 


Figure  6.  Assessment  of  the  noise-based  estimation  of  time  transfer  uncertainty  for  areal  data  sequence.  Cross  (+)  represent 
simulated  time  offsets  (eg.  over  5  days).  The  boxes  and  its  error  bar  represent  the  daily  interpolated  time  offset  and  its  standard 
deviation.  D ashed  lines  represent  the  noi se-base d  standard  devi ati on  around  the  interpol ation  solid  line 
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Questions  and  Answers 

JOHN  DAVIS  (National  Physical  Laboratory,  UK):  Did  you  do  any  statistics  just  looking  at  GPS 
measurements  at  the  single  epoch?  That  is,  measurements  that  were  all  taken  at  the  same  time,  presuming 
that  you  were  using  a  multi-channel  GPS  receiver?  Are  you  using  a  multi-channel  GPS  receiver  to  make 
your  measurements? 

MAHMOUD  ADDOUCHE:  I  use  the  BIPM  schedule.  So  I  track  one  satellite  at  one  time. 

DAVIS:  The  point  was  that  if  you  are  using  a  multi-channel  receiver  and  you  are  making  your 
measurements  at  the  same  epoch,  then  if  you  look  at  the  statistics  on  that,  you  will  be  independent  of  the 
clock  noise;  whereas,  if  you  are  making  measurements  at  successive  epochs  with  a  single-channel 
receiver,  you  are  going  to  get  a  mixture  of  clock  and  time  transfer  noise.  This  is  an  easy  way  of 
distinguishing  the  two,  which  might  help  you. 

ADDOUCHE:  It  will  be  interesting,  but  I  didn’t  see  this. 

DAVIS:  There  should  be  plenty  of  data  available  for  you  to  have  a  look  at.  We  can  provide  you  with  that. 
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